#include <iostream>
#include <new>
#include <math.h>
#include <stdlib.h>

const double MASS = 9.10956e-31;  /* Kg */
const double Q = 1.6e-19;	     /* Coulomb */
const double PLANK = 6.6262e-34;   /* Js */
const double PLANKUT = 1.05459e-34;  /* Js */
const double PI	= 3.14159;
const double KB	= 1.38062e-23;  /* J over K */
const double VT	= 0.024;	     /* eV */
const double BOHRAD = 5.29177e-11;  /* meters */
const double EPS0 = 8.854e-14;    /* Farad over cm */

const double N2DEG = (MASS/PLANKUT)*(VT*Q/(PI*PLANKUT));
const double N3DEG = 1e-6*(VT*Q*sqrt(Q)/PLANKUT)*(MASS/PLANKUT)*(sqrt(MASS)/(PI*PI*PLANKUT));

const int MAX_X = 500;
const int MAX_Y = 250;
const int MAX_T = 50;
const double L = 1.0e-4;

const double KKK = 26.19168482;	/* 10^-18*2*q*me  h^2 in nm */
const double LNTHCONV =	1.0e7;	/* Factor for conversion from cm to nm */

/*Max number of quantum steps allowed */
/*A quantum step is a discontinuity in the /* 
/* energy band that can be part of a quantum well */
const int MAX_TRANSITIONS = 10;

